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Abstract 

Results on soft and hard diffraction are briefly reviewed and placed in a QCD per- 
spective using a parton model approach. Issues addressed include factorization, scaling 
properties, universality of rapidity gap formation, and unitarity. Predictions for differen- 
tial cross sections of processes with multiple rapidity gaps are presented with examples 
for the Tevatron and the Large Hadron Collider. This paper is an expanded version of 
a paper delivered at "Snowmass-2001" 



1 Introduction 

Experiments at p(p)-p and e-p colliders have reported and characterized a class of events 
incorporating a hard (high transverse momentum) partonic scattering while carrying 
the characteristic signature of diffraction, namely a leading (anti)proton and/or a large 
"rapidity gap" , defined as a region of pseudorapidity, i] = — ln(tan|), devoid of particles 
(see Figs. 1 and 2). The rapidity gap is a non-perturbative phenomenon presumed to be 
caused by the exchange of a Pomeron , whose generic QCD definition is a color-singlet 
combination of quarks and/or gluons carrying the quantum numbers of the vacuum. 
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FIG. 1: Dijet production in pp single (a) and double (b) FIG. 2: Diagram for diffractive deep inelastic 
diffraction, and in double Pomeron exchange (c). scattering at HERA. 



The interplay between soft and hard processes in hard diffraction is of particular the- 
oretical importance due to its potential for elucidating the transition from perturbative 
to nonperturbative QCD. Phenomenological models proposed for hard diffraction have 
been only partially successful in describing the data. A QCD-based theoretical descrip- 
tion is still not available. This is not surprising, since diffraction invariably involves 
non-perturbative effects associated with the formation of rapidity gaps. The theoretical 
community has thus far paid very little attention to soft diffractive processes, i.e. those 
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which do not have a hard partonic scattering. Yet, experiment has shown that soft pro- 
cesses have "a lot to say" about rapidity gap formation, which is key to understanding 
diffraction. In this paper, we examine the regularities observed in both soft and hard 
diffractive processes using a parton model approach, which places diffraction in a QCD 
perspective. 



2 Rapidity gaps 

The exchange of a gluon or a quark between two hadrons at high energies leads to events 
in which, in addition to whatever hard scattering may have occurred, the entire rapidity 
space is filled with soft (low transverse momentum) particles (underlying event). The 
soft particle distribution is approximately flat in (pseudo) rapidity. The flat dN/dr] shape 
is the result of x-scaling |J of the parton distribution functions of the incoming hadrons. 
Rapidity gaps can be formed in any inelastic non- diffractive (ND) event by multiplicity 
fluctuations. Assuming Poisson statistics, the probability for a gap of width Arj within 
a ND event sample is given by 

P™(&v)=pe- pAr> (1) 

where p is the average particle density per unit rj (the total probability for any gap 
is normalized to unity). As seen in ([!]), the probability for rapidity gaps formed my 
multiplicity fluctuations is exponentially suppressed. In contrast, since no radiation 
is emitted by the acceleration of vacuum quantum numbers, Pomeron exchange leads 
to large rapidity gaps whose probability is not exponentially dumped. Therefore, large 
rapidity gaps are an unmistakable diffractive signature and can be considered the generic 
definition of diffraction. 



3 Soft diffraction 



Soft diffraction has been traditionally treated theoretically in the framework of Regge 
theory. For large rapidity gaps (Ar/ > 3), the cross sections for p(p)-p single and double 
(central) diffraction can be written as M 
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where a(t) = 1+e+a't is the Pomeron trajectory, (3{t) the coupling of the Pomeron to the 
(anti)proton, and k = g(t)//3(0) the ratio of the triple- Pomeron to the Pomeron-proton 
couplings. The above two equations, which are based on factorization, are remarkably 
similar. In each case, there are two factors: 
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the first, which is 



the second, which has the form 
squared s'. 



{exp[(e + a't)Arj}} 2 and thus depends on the rapidity gap, 

(s'Y of a total cross section at c.m.s. energy 
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In single diffraction (SD), the second factor is interpreted as the Pomeron-nucleon to- 
tal cross section, while the first factor as the square of the elastic scattering amplitude 
between the diffractively excited nucleon state and the other nucleon. The similarity 
between the SD and double diffraction (DD) equations suggests an interpretation of 
the first factor in DD as the elastic scattering between the two diffractively dissociated 
nucleon states. However, the interpretation of the second factor in DD is not as straight- 
forward as in SD, but does acquire a similar meaning in the parton model, as discussed 
in the next section. Since the second factors in (0) represent properly normalized cross 
sections, the first factors may be thought of as the rapidity gap probability and renor- 
malized to unity. A model based on such a renormalization procedure |5|, U has yielded 
predictions in excellent agreement with measured SD and DD cross sections, as seen in 
Figs. 3 and 4. 
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FIG. 3: The pp total SD cross section exhibits an s- 
dependence consistent with the renormalization proce- 
dure of Ref . [5] , contrary to the s 2e behaviour expected 
from Regge theory (figure from Ref. [5]). 
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FIG. 4: The pp total DD (central gap) cross section 
agrees with the prediction of the renormalized rapidity 
gap model [6], contrary to the s 2e expectation from 
Regge theory (figure from Ref. [4]). 



4 Factorization, scaling and unitarity 

The renormalized rapidity gap probability is by definition energy independent and thus 
represents a scaling behaviour, which is most prominently displayed in the form of the 
SD differential cross section in terms of the square of the mass of the diffractively excited 
nucleon state, M 2 . The variable M 2 is related to Ar) and to the fractional momentum 
loss of the leading nucleon: M 2 = s£ = se _Ar? . The unrenormalized SD cross section at 
t = has the form ^ffi\ t=0 oc ( M s 2 y +e ■ This form leads to a total diffractive cross section 
~ s 2e , which grows faster than the total nucleon-nucleon cross section (~ s e ) and would 
exceed it at ^/s «2 TeV violating unitarity (see Figs. 3 and 4). The renormalization 
procedure leads to the energy independent form ^fM-\ t=0 cx tmT\t+i , which avoids the 
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violation of unitarity. This form is brilliantly supported by the data, as can be seen in 
Fig. 5. Thus, it appears that nature preserves unitarity by favoring "M 2 -scaling" over 
factorization. 
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Fig. 5: Cross sections d 2 asD/dM 2 dt for p + p(p) — > p(p) + X at t = and ^/s = 14, 
20, 546 and 1800 GeV, multiplied by [/^pS 2 ^^ • 9ipipip) /(16tt * V(s))]" 1 , where 
N(s) is the integral of the rapidity gap probability, are compared with the renormalized 
gap prediction of l/(M 2 ) 1+e using e = 0.104 0. The dashed curves show the standard 
Regge-theory predictions. The t = data were obtained by extrapolation from their 
t = —0.05 GeV 2 values after subtracting the pion exchange contribution M. 



5 A parton model approach 

In the parton model, the pp total cross section is basically proportional to the number 
of partons in the proton, integrated down to x — s J s, where s Q is the energy scale for 
soft physics. The latter is of C((M T ) 2 ), where (M T ) ~ 1 GeV is the average transverse 
mass of the particles in the final state. Expressing the parton density as a power law in 
1/x, which is an appropriate parameterization for the small x region responsible for the 
cross section rise at high energies, we obtain 

f 1 dx ( s 

J (so A) x 1+n \s 

The parameter n is identified as the e = a(0) — 1 of the Pomeron trajectory 0. 
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According to the optical theorem, the form s e of the total cross section is that of the 
imaginary part of the forward (t = 0) elastic scattering amplitude. Note that In s is the 
(pseudo)rapidity region Arj in which there is particle production. The full parton model 
amplitude can be written as 

lmi(t,Ay)^ e ^ + a '^ A y (4) 

where we have added to e the term oit as a parameterization of the ^-dependence of the 
amplitude. 

One may now understand expressions (Q) for the SD and DD cross sections as follows: 
the second term in brackets represents the nucleon-nucleon total cross section at a sub- 
energy squared s', multiplied by a factor k, which may be interpreted as "the price 
to pay" for the color matching required to enable the formation of a rapidity gap, 
while the first term represents the amplitude squared of the elastic scattering between 
the diffractively dissociated nucleon and the other nucleon in SD, or between the two 
diffr actively dissociated nucleons in DD. Since the cross section represented by the second 
term is properly normalized, the first term simply represents a rapidity gap probability 
distribution and should be normalized to unity. This result is equivalent to renormalizing 
the rapidity gap probability of the Regge model, as proposed in ||, || . 



6 Multiple rapidity gaps in diffraction 

The parton model approach used above to calculate the SD and DD differential cross 
sections lends to easy generalization to events with multiple rapidity gaps. Here, vie 
outline the procedure for calculating the differential cross section for a 4-gap event (see 
Fig. 6). 
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Fig. 6: Topology of a 4-gap event in pseudorapidity space. 

The calculation of the differential cross section is based on the parton-model scatter- 
ing amplitude of Eq. (|4|). For the rapidity regions A?]-, where there is particle production, 

the t = parton model amplitude is used and the sub-energy cross section is ~ e ^ . 
For rapidity gaps Ar) which can be considered as resulting from elastic scattering be- 
tween diffractively excited states, the square of the full parton-model amplitude is used, 
g 2(e + a ti)Ar)i^ anc j form factor (3 2 (t) is included for a surviving nucleon. The gap 
probability (product of all rapidity gap terms) is then normalized to unity, and a color 
matching factor k is included for each gap. 
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For the 4-gap example of Fig. 6, which has 10 independent variables Vi (shown below 
the figure), we have: 



= Pgap x <x(sub - energy) 

(3 2 (0) ■ e eA y' 



a (sub — energy) = k 



(Ay' = E-=i Mi, 



• P gap = N gap x n| =1 [e^)^] 2 x \P(h)P(U)] 2 = N gap • e 2e M • fW\£i (M = 

Eii M) 

where N gap is the factor that normalizes P gap over all phase space to unity. 

The last equation shows that the renormalization factor of the gap probability depends 
only on s (since Ai] max = Ins) and not on the number of diffractive gaps. Thus, the 
ratio of two-gap to one-gap cross sections is expected to be ~ k, with no additional 
energy dependent suppression for the second gap. Below, we discuss two specific two- 
gap processes which can be studied at the Tevatron. Studies of processes with more 
than two gaps will have to await the commissioning of the LHC 

6.1 Double Pomeron exchange (DPE) 

The double Pomeron exchange process is shown below for pp collisions (we use rapidity 
y and pseudorapidity r\ interchangeably): 

Ayi ^ ^ Ay 2 

y c 

t\ Ay = Ay 1 + Ay 2 t 2 

Fig. 7: p + p -> p + GAP + X + GAP + p 
The Regge theory predictions for SD and DPE are: 
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Following the procedure outlined above, the DPE differential cross section can be 
written as 
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The renormalization factor N^pe is obtained by integrating the gap probability P gap 
over ti, t 2 , Ay, y c . The limits of integration are: 
< -ti < oo; < -t 2 < oo 

2.3 < Ay < m(s/s ), where sq — 1 GeV 2 and Ay = 2.3 corresponds to £ = 0.1 
-\{Ay - 2.3) < y c < ±(Ay - 2.3) 

For SD, the gap probability is the first term in the SD equation and N$d is obtained by 
integrating it over Ay within the region 2.3 < Ay < ln(s/so). 

For numerical evaluations we use § 1 a(t) = 1.104 + 0.25t, (3(0) = 4.1 mb* (6.57 
GeV -1 ) and k = 0.17. The ratio of DPE to SD rates in pp collisions with a leading 
final state antiproton of fractional momentum loss £p = 0.065, the average value of the 



inclusive p-triggered data samples in the CDF publications of Refs. |TIJ [TT]] , and for 
£ p > 0.02 is predicted to be 



jdDPE 



fp -=o.o65, fp >o.o2 = 0.21 (0.17) ± 10% at 1800 (630) GeV 



where the error is due the the uncertainty on the factor k. This ratio does not depend 
strongly on £(p) within the range 0.035 < £(p) < 0.95 of the CDF event samples and 
therefore represents a realistic prediction. 



6.2 Single diffraction with a second gap (SDD) 

An interesting two-gap process is pp single diffraction with a leading final state antiproton 
and a rapidity gap within the rapidity space allocated to the diffraction dissociation 
products of the proton (see Fig. 7). 
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FIG. 8: p + p -> p + GAPi + X + GAP 2 + Y 

The second gap, GAP 2 , may be thought of as due to W-p double diffraction. Thus, 
this process is a combination of single plus double diffraction and hence we represent it 
by "SDD". 
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To calculate the SDD cross section we start with the Regge predictions for SD, DD 
and SDD: 
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Following the rules of the parton model approach, the SDD cross section may be written 
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where In s" = In s — Ayi — A?/2 



Changing variables, Ayi, Ay 2 =>■ Ay = Aj/x + Ay 2 , J/rf, where yd = \Ay\ is the center 
of the single diffractive cluster, the gap probability becomes 



Pga P (t 1 ,t 2 , Ay, y d , y c ) = 2k 



16tt 
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where -Fi(ii) is the antiproton form factor ||. Integrating P gap over ti,t 2 ,y c ,yd, Ay 
yields the renormalization factor N$dd- The limits of integration are: < —t\ < oo; 



-(Ay~2y d ) 



< -t 2 < e A ^; -\n^s + 2y d + \{Ay - 2y d ) < y c < In y/S - \{Ay - 2y d 



0<y d < \Ay- 2.3 < Ay < In (1-5 a ^L a , where s = 1 GeV 2 . 

Experimentally, in order to be able to detect a large rapidity gap within a single- 
diffractive cluster of particles, the cluster must extend over a large part of the rapidity 
space covered by the detector. Therefore, we evaluate the SDD to SD fraction for 
events with a relatively large £(p)-value. For £ p = 0.075, which for yfs =1800 (630) 
GeV corresponds to sub-energy yfs' m 500 (170) GeV, the predicted SDD fraction for 
Ay 2 > 3 is 



•ptwo—gap 
-^one—gap 



€p -=o.o75,a s/2 >3 = 0.26 (0.20) ± 10% at 1800 (630) GeV 



where the error is due the the uncertainty on k. This prediction can be compared with 



experiment using the existing CDF inclusive SD data of Refs. |L0], |TT . 
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7 Hard diffraction 



The central issue in hard diffraction is the question of QCD factorization: can hard 
diffractive cross sections be obtained as a convolution of "diffractive structure functions" 
(DSF) with parton-level cross sections? This question was addressed decisively by a 
comparison |T(| between the DSF measured by CDF in dijet production at the Tevatron 
and the prediction based on parton densities extracted from diffractive DIS at HERA. 
The DSF at the Tevatron was found to be suppressed relative to the prediction from 
HERA by a factor of ~ 10 (see Fig. 9). This result confirmed previous CDF results 
from diffractive W dijet |L3| and 6-quark |14| production at ^=1800 GeV, and 
was corroborated by more recent CDF results on diffractive J/ip [ITS] production at 1800 



GeV and dijet production at 630 GeV [ I I 




Fig. 9: The diffractive structure function Fjj(f3,Q 2 ), where j3 is the momentum 
fraction of the parton in the Pomeron, extracted from CDF diffractive dijet production 

data in pp collisions at y/s =1800 GeV at the Fermilab Tevatron, is compared with 
expectations from parton densities extracted from diffractive deep inelastic scattering 
by the HI Collaboration at the DESY ep collider HERA (figure from Ref. [0). 

Although factorization breaks down severely between HERA and the Tevatron, it nev- 
ertheless holds within the HERA data and within the single-diffractive data at the 
Tevatron at the same center of mass collision energy. This is demonstrated by the fact 
that the gluon parton distribution function (PDF) derived from DIS adequately de- 
scribes diffractive dijet production at HERA |T6| , while at the Tevatron a consistent 
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gluon PDF is obtained from the measured rates of diffractive W, dijet, fe-quark and J/ip 
production |TI| . 

Factorization was also tested at the Tevatron between the structure functions mea- 
sured in single-diffractive and double-Pomeron exchange (DPE) dijet production at 



y/s =1800 GeV fll7fl . The ratio of the DPE to SD structure functions was found to 
be larger than that of the SD to ND ones by a factor of 5.3 ± 2.0. This result represents 
a breakdown of factorization within Tevatron data. However, we note that DPE is a 
two-gap process and, according to the parton model approach to diffraction that we 
have presented, should not be suppressed (except for kinematical edge effects) relative 
to SD. The above result, within the experimental uncertainty, confirms the prediction 
of the parton model approach. 

Finally, CDF has reported that the f3 and £ dependence of the diffractive structure 
function at y/s = 1800 GeV in the region 0.035 < £ < 0.095 obeys /3-£ factorization for 
(3 < 0.5 The observed ~ £ _1 dependence shows that Pomeron-like behaviour extends to 
moderately high £ values in diffractive dijet production. Such behaviour is expected in 
the parton model approach for hard diffraction, where the Pomeron emerges from the 
quark-gluon combination of two partonic exchanges, one on a hard scale that 

produces the dijet system and the other on a soft scale that neutralizes the color flow 
and forms the rapidity gap [ |T8|j . 



8 Conclusion 

Soft hadronic processes have traditionally been treated theoretically in the framework of 
Regge theory. By postulating a Pomeron with a trajectory of intercept a(0) > 1, Regge 
theory correctly predicts several features of soft physics including the rise of total cross 
sections with increasing energy, the shrinking of the forward elastic scattering peak, 
elastic to total cross section ratios, and the shape of the single diffraction differential 
cross section. However, the predicted ~ s 2e rise of the single diffractive cross section 
is faster than the ~ s e rise of <jt, so that if <7sd is normalized to the experimental 
value at ^/s ~ 22 GeV the Regge prediction would exceed <jt at y/s rs 2 TeV violat- 
ing unitarity [[J. The violation of unitarity is averted by renormalizing the "Pomeron 
flux" H to unity, which is equivalent to renormalizing the rapidity gap probability of 
Eq. @. This prescription leads to excellent agreement with data. It is interesting that 
the renormalization procedure, which violates Regge factorization, results in an energy 
independent da,sD/dM 2 cross section, which represents a scaling behaviour; moreover, 
it can also be applied to double-diffraction ||. In the present paper, these results are 
obtained in a parton model approach, which has the added advantage that it can be 
generalized to describe processes with multiple rapidity gaps H. The most interesting 
conclusion from this approach is that in multigap events the renormalization factor de- 
pends only on the s-value of the collision, so that the same suppression factor, otherwise 
known as "gap survival probability", applies to a process regardless of the number of 
gaps it contains. The ratio of rates of two-gap to one-gap production in a diffractive 
process is predicted to be equal to the ratio of the PPP to P-p couplings, neglecting 
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phase space edge effects. Thus, this factor is interpreted as being the "price to pay" 
for the color-matching required to form a color-singlet with vacuum quantum numbers. 
The parton model approach can also be applied to hard diffraction (not discussed here 
for lack of space) and explains the breakdown of diffractive QCD factorization between 
HERA and the Tevatron [If . 
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